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ABSTRACT 

We study global non-axisymmetric oscillation modes trapped near the inner boundary 
of an accretion disc. Observations indicate that some of the quasi-periodic oscillations 
(QPOs) observed in the luminosities of accreting compact objects (neutron stars, black 
holes and white dwarfs) are produced in the inner-most regions of accretion discs or 
boundary layers. Two simple models are considered in this paper: The magnetosphere- 
disc model consists of a thin Keplerian disc in contact with a uniformly rotating 
magnetosphere with and low plasma density, while the star-disc model involves a Ke- 
plerian disc terminated at the stellar atomosphere with high density and small density 
scale height. We find that the interface modes at the magnetosphere-disc boundary 
are generally unstable due to Ray leigh- Taylor and/or Kelvin-Helmholtz instabilities. 
However, differential rotation of the disc tends to suppress Rayleigh- Taylor instability 
and a sufficiently high disc sound speed (or temperature) is needed to overcome this 
suppression and to attain net mode growth. On the other hand, Kelvin-Helmholtz 
instability may be active at low disc sound speeds. We also find that the interface 
modes trapped at the boundary between a thin disc and an unmagnetized star do not 
suffer Rayleigh- Taylor or Kelvin-Helmholtz instability, but can become unstable due 
to wave leakage to large disc radii and, for sufficiently steep disc density distributions, 
due to wave absorption at the corotation resonance in the disc. The non-axisymmetric 
interface modes studied in this paper may be relevant to the high-frequency QPOs 
observed in some X-ray binaries and in cataclysmic variables. 

1 INTRODUCTION 

Quasi-periodic variabilities have been observed in the timing data of various types of accreting objects. Several types of 
quasi-periodic oscillations (QPOs) are observed in X-ray binaries with accreting black holes (BHs) or neutron stars (NSs) 
(e.g., Remillard & McClintock 2006; Van der Klis 2006). Oscillations are also seen in the outbursts of accreting white dwarf 
(WD) systems (e.g., Patterson 1981; see Warner 2004 for a review). 

In accreting NS and BH X-ray binaries the observed QPO frequencies (40 — 450 Hz for the high-frequency QPOs in the 
BH systems and > 300 Hz for kHz QPOs in the NS systems) imply a source close to the central compact object where the 
Keplerian orbital frequencies are high. Since the BH systems lack a hard surface where oscillations may occur, it is likely that 
the source of the variability is in the inner regions of the disc itself or in some interface regions between the disc and the 
plunging flow. Gilfanov et al (2003), however, found that, based on spectral analysis of the disc emission components, the 
quasi-periodic variability in Low Mass NS X-ray Binary systems are most likely caused by variations in the disc boundary 
layer, rather than the disc itself. 

In Cataclysmic Variables (CVs) the Dwarf Nova Oscillations (DNOs) seen during during outbursts have frequencies 
roughly corresponding to the Keplerian rotation rate at the WD surface (e.g., Patterson, 1981; Warner, 2004; Knigge et al, 
1998), which imply an origin at or near the inner disk boundary. 

Several models involving accretion disk boundary dynamics have been proposed in different contexts. Popham (1999) 
studied the effect of a non-axysymmetric bulge at the optically thick to optically thin transition radius as a model for DNOs. 
Piro & Bildsten (2004) examined the surface wave oscillations that would occur within the thin equatorial belt around a 
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non-magnetized WD formed by the accretion spreading layer, while Warner & Woudt (2002) considered accretion onto a 
slipping belt. In the context of accreting magnetic (neutron) stars, Arons & Lea (1976) and Eisner & Lamb (1977) considered 
the interchange instability at the magnetosphere boundary. Spruit & Taam (1990) and Spruit, Stehle & Papaloizou (1995) 
investigated the stability of thin rotating magnetized discs. Of particular relevance to the present paper is the work of Li & 
Narayan (2004), who examined a simplified cylindrical model of the Rayleigh- Taylor and Kelvin- Helmholtz instabilities at 
the boundary between a magnetosphere and an incompressible rotating flow. There have also been a number of numerical 
simulations of the interface at the magnetosphere-disc boundary (see Romanova et al., 2008 and Kulkarni & Romanova, 2008 
and references therein). 

In this paper we study global non-axisymmetric oscillation modes confined near inner boundary of the accretion disc 
(interface modes). We consider two simple models. The first model involves the magnetosphere-disc boundary similar to the 
model of Li & Narayan (2004) : we consider an uniformly rotating incompressible magnetosphere with low gas density (where 
magnetic pressure dominates), which truncates a thin barotropic accretion disc (where gas pressure dominates). This situation 
may arise from magnetic field build up due to accretion (e.g., Bisnovatyi-Kogan & Ruzmaikin, 1974, 1976; Igumenshchev et 
al., 2003; Rothstein & Lovelace, 2008) or by the magnetosphere of a central (neutron) star. Unlike Li & Narayan (2004), who 
restricted their model to incompressible fluid, our discs are compressible and we show that because of the differential rotation 
of the disc, finite disc sound speed plays an important role in the development of the instability of the interface modes. 

In our second model we examine the interface modes for accretion onto a non-magnetic stellar surface. Though the 
structure of the boundary layer is non-trivial and may affect boundary modes (see, e.g., Carroll et al., 1985; Collins et al., 
2000), we consider the instabilities for a thin disc truncated by a sharp transition to a dense uniformly rotating stellar 
atmosphere. This simplified model may provide insight for modes with characteristic radial length scale much greater than 
the radial length scale of the boundary layer. 

In Section 2 we describe the basic setup for the magnetospheric boundary model, and in Section 3 we discuss the resulting 
interface mode instabilities. We describe the star-disc boundary and analyze its possible instabilities in Section 4. We then 
conclude in Section 5 with a discussion of possible applications of our findings. 



2 MAGNETOSPHERE-DISC SETUP 

We begin by considering a simplified model of the magnetosphere-disc boundary similar to the one considered by Li & Narayan 
(2004). The magnetic field is assumed to be negligible in the disc region (r > ri n ), while the magnetosphere region (r < ri n ) 
is assumed to be incompressible and have low density compared to the disc region, with purely vertical magnetic field. Unlike 
Li & Narayan (2004), who assumed infinite sound speed in the disc, our disc has sound speed c s much less than the disk 
rotation speed rQ. 

In terms of the vertically integrated density (E), pressure (P), magnetic field (B) and fluid velocity (u) the ideal MHD 
equations are: 

5g+V.(Eu) = (1) 
^ + (u-V)u = -Ivn-V$ + iT (2) 



dB 
~dt 



= V x (u x B). (3) 



where II = P + B 2 /8n is the total pressure, T = ^L(B ■ V)B is the magnetic tension, and $ is the gravitational potential due 
to the central object (e.g. Fu & Lai, 2008). Using cylindrical coordinates (r, tj>, z), we consider the case where the magnetic 
field is purely poloidal and B = B z z in the disc plane, which gives T = 0. We assume an axisymmetric background flow with 
fluid velocity u = rQ(r)<f>. The unperturbed flow satisfies the condition 

1 dn d$ 

^-E* = dr-" V (4) 
The linearized equations of |T} and ([2| with perturbations of the form g™*'- 1 "' (assuming no vertical dependence) take 
the form: 

- iuST + -^-(Tr8u r ) + ^^-Sua = 0, (5) 
r or r 

—icoSur — 2i}Su,p = — Pcff^P — ^r-^-SU , (6) 

L L or 

2 ^ m 

-iuSus + —Su r = - — 511 , (7) 
III T,r 

where Cu — uj — mQ. is the Doppler shifted frequency, « = [^^(r 2 ^)] is the radial epicyclic frequency, and ST,, 5T1 and <5u 
are the Eulerian perturbations of the fluid variables. Additionally, assuming a barotropic flow we have 
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Figure 1. The wavefunctions for the interface mode for m = 4, c s = O.lrQ, = 1/99 and Q^/Q- ln = 1, with the mode frequency 

u>/Q- ln = 4.275 + 0.1914i, where Q ln = f2(ri n +) ~ Ofe(r"i n ). The real components are shown in solid lines, while the imaginary components 
are dashed lines. Note that for the interface modes only the outer Lindblad resonance (rpLi denoted by the dotted line) exists outside 



5 E=V=>I-i-B.*B). (8) 



with the sound speed c s = (dP/dS) 



1/2 



2.1 The Magnetosphere 



In the inner, magnetically dominated region (r < n n ), we assume the flow to be incompressible, and have uniform rotation 
(fl = Q- = const) and uniform surface density (E = E_ = const). Equations © - (0 then reduce to 



ld_ 

r dr 



(ErSur 



imTi 



Suj, = 



r 



i_d_ 

EcV 
im 
~Er 



5n 

SIL 



As in Li & Narayan (2004) we define W = rSu r and find 8u^ = (i/m)dW/dr and 



ld_ f dW\ n? 



r dr 



dr I 



1 



r d 
muj dr 



K 

20 



W = 



(9) 
(10) 

(11) 
(12) 



For uniform rotation, k = 20, equation (|12f) has the solution W oc r ±m . Since r < r ln , we take the positive sign to be the 
physical solution so that the perturbation falls off away from the interface. Thus the exact solution for the r < ri n region is 

m — 1 



8u r = Su r (rin) 



(13) 



2.2 The Disc 

In the disc (r > n n ), we take the magnetic field to be small, such that P ^> B 2 /(8ty), and the angular velocity of the 
unperturbed flow to be nearly Keplerian, such that f2(r) « Qk(r) = \J \^§r- Rewriting equations (J5]l - ([7]) we have 

.. E 1 d -_. . . imE . 

- iuj—5h+ - — (ErSur) H Su^ = 0, (14) 

cj r or r 
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-iu)8u r 



2Q,8u ( j, 

2 



dr 
im 
r 



5h , (15) 
5h , (16) 

E E < 17 > 
is the enthalpy perturbation. Eliminating the velocity perturbations in favor of the enthalpy, we obtain the second order ODE 
for the enthalpy perturbation in the disc, 



where 

2 5E SP 
on = c, — — = — — 



dr 2 



d_ 

dr 



In- 



D \ d 



rE / dr 



d , fiE 
-r- m "TT- 



m 



D 
7? 



<5ft = 



(18) 



where D = k 2 — uj 2 . For concreteness we will assume a power-law disc surface density profile Ear p . 



2.3 Matching Conditions Across the Interface 



The matching conditions across the interface at r ln between the magnetosphere and the disc region are given by demanding the 
continuity of the Lagrangian displacement in the radial direction £ r = i5u r /u, and the total Lagrangian pressure perturbation 
All = (511 + £ r across the boundary. The former gives 



i8u r -\ 



iSu r - 



(19) 

where the subscript "±" implifes that the quantities are evaluated at r = n n ±. The total Lagrangian pressure perturbation 
for r = ri n _ is given by 

2 



An, 



= E. 



2m£l 
2rflQ 



+ 



i dfl 



ru) dr 

9cs- + 





iruj 


w + — 




m A 
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Hi 



(20) 



In the disc region, we have 
An+ = AP+ = E 



i&u r 1 dP 
5h+—- — 

uj E dr 



= E. 



/ u>8h 
\iSu r 



PCs 

r 



iSu r 
UJ+ 



The condition AII4 



/ CjSh pc 2 
\iSu r 



AII_ then gives 
(2rQ.Hi 



= E_ 



\ m 



9eS 



rui 
m 



(21) 



(22) 



3 INTERFACE MODES AT THE MAGNETOSPHERE-DISC BOUNDARY 

Perturbations mainly confined to the magnetosphere-disc interface can become unstable due to Ray leigh- Taylor or Kelvin 
Helmholtz instability. In order to calculate the growth rates, we must solve the eigenvalue problem given by equation (|18[) 
with an outgoing wave boundary condition at some outer radius, and equation (|22p at the interface radius n n . 



3.1 Numerical Solution 

We adopt the radiative outer boundary condition in the outer wave zone of the disc, such that far from the outer Lindblad 
resonance radius tql (where to — mQ = k) we have the solution of the form: 



ShxAexp^iJ kdr J , (23) 

with A = (L>/rEfc) 1/2 and k = (~D/c 2 ) 1/2 (see Tsang & Lai 2008, Lai & Tsang, 2008). This gives the boundary condition at 
r = r ou t > roh- 
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Figure 2. Real and imaginary frequencies for interface modes for various m as a function of the disc sound speed c s . The left panels 
are for S_/S+ = 0, while the right panels are for E_/E_|_ = 1/9, both with Q—/Q ln = 1. The solid lines show the eigenfrequencies for 
m = 5 modes, the dash-dotted lines for m = 4, the long-dashed lines for m = 3, the short-dashed lines for m = 2, and the dotted lines 
for m = 1. 




Sh'(r out ) = Sh(r out ) Jjfc+I^l) . (24) 

V ' font 

We adopt (|22p as the inner boundary condition for the disc and solve the eigenvalue problem using a standard shooting 
method (Press et al 1998). For the numerical solutions below, the density profile of the disk was assumed to be E ~ r ~ 3 / 2 S o 
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that corotation absorption plays no role in determining the mode stability (Tsang & Lai 2008). An example wavefunction for 
an interface mode is shown in Figure 1, for typical disc parameters. 

The numerical eigenvalues are shown in Figure 2 and Figure 3 for various disc and magnetosphere parameters, for 
m — 1,2, ... ,5. 



3.2 Discussion of Numerical Results 

Figure 2 shows the complex eigenvalues as a function of sound speed c s , for density contrasts corresponding to E_ = and 
E_ = |E+, with magnetosphere rotation rate equal to the Kepler frequency at the interface = 0(rin+) ~ £2fc(Tm)]. For 
this case we see that there exists a cutoff in the disc sound speed below which no growing interface modes are found. This 
arises from the stabilizing effect of the background differential rotation, and can be understood as follows. 
Setting E_ = and rewriting (|22[) in terms of the radial velocity perturbation Su r , we have: 

£> 2 d / ru 2 m 



where g c a+ = pc 2 /r [see Eq. (J4j>] and £ = k 2 /(2J7) is the vorticity, and where all quantities are evaluated at the interface 
n n +. For the wave frequencies of interest, the waves are evanescent in the region of the disk just outside the interface. Let 
k = —8u' r /Su r > 0. Equation ([25} can be solved in terms of k, giving 



1 g cff+ m 2 


m 2 C, 2 


r7 


4 7 2 



where y = kr — 1. The terms inside the square root correspond to the mode growth due to Ray leigh- Taylor instability and 
the suppression due to vorticity, respectively. With g B s+ = pc 2 /r we find the critical sound speed 



C 2 r 2 

**«V^"' (27) 

above which the perturbations will be unstable. 

Figure 3 shows cases where the inner region is uniformly rotating at an angular frequency of one half the Kepler frequency 
at the interface [fi_ = 0.5f2(ri n +)]- When E_ > this leads to the development of the Kelvin-Helmholtz instability, and both 
this and the Rayleigh- Taylor instability play a role in the mode growth. In the Appendix we derive the expression for 
the plane-parallel Rayleigh- Taylor and Kelvin-Helmholtz instabilities for a compressible upper region (with density p+ and 
horizontal velocity u+), and incompressible lower region (with density p_ and horizontal velocity For p_ <C p+ we have 
u) w ku+ ± iuJi where k is the horizontal wavenumber and 



kk(u+ -u-) 2 +gk= ^ H + W 2 T . (28) 

Here k « [y/1 + H 2 k 2 — f] , g is the acceleration due to gravity in the vertical direction, and H z is the vertical scale 
height in the upper region. The Kelvin-Helmholtz term is approximately 

WKH~tt(«+-tt-f^j. (29) 

For kH z 3> 1 this reduces to the incompressible limit with o^h ~ k 2 (u+ — u~) 2 p_ / p+. For kH z < 1 we have w| H w 
(H z k)k 2 (u+ — U-) 2 p- / p+, a factor of H z k smaller than the incompressible result. 

For the rotating system under consideration, the imaginary part of the mode frequency can be written schematically as 
[cf. equation ([28}] 



+ <4 T +o4 rt . (30) 



We also have H z ~ c 2 /g c ff ~ r and k ~ m/r so kH z ~ m. Thus wkh ~ 

(Afi) 2 E_/E+, and ojkh depends weakly on 
sound speed. On the other hand, from equation (I26|l we see that the vorticity suppresses mode growth through the term 



vort = — m ^ /(4r/ ). For sufficiently small c s , equation (|18[) indicates 8h oc e _Kr ' Cs , i.e. A: ~ k/c s . Therefore the vorticity 
term scales with sound speed as Lo 2 OIt ~ —m 2 c 2 /r 2 , and can be dominated by the Kelvin-Helmholtz term for small enough 
sound speed. In the left panel of Figure 3, the mode growth (u>i > 0) for small c s is mainly driven by the Kelvin-Helmholtz 
instability. For m ^ 5 the sound speed ranges where w^t an d ^kh dominate over w 2 OTt overlap, and hence the critical sound 
speed in equation (|27p is not relevant. For larger values of E_ these regions can overlap for all m. 
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Figure 4. Real and imaginary frequencies for the interface modes of various m (denoted as in Figure 2) as a function of sound speed c s 
with ri n = nsco f° r the pseudo-newtonian GR potential. There is no sound speed cutoff for unstable modes as the vorticity is zero at the 
interface. The left panels show the eigenvalues for S_ = and f2_ = Q(r; n _|_) and the right panel shows the eigenvalues for S_ = ijg£+ 
and f2_ = 0.5f2i n 



3.3 Effect of a Relativistic Potential 



While in Section 3.2 and other sections of the paper we focus on Newtonian discs, it is of interest to consider how general 
relativity may modify our results. The effect of general relativity can be approximated by using the pseudo-Newtonian 
Paczynski & Wiita (1980) potential: 



GM 
r — rs' 



(31) 



with rs = 2GM/c? the Schwarzschild radius. This gives the Keplerian orbital frequency (flk) and epicyclic frequency (/t) as: 



I— 

r dr 



1/2 



GM 



rs 



— (r fife) 

r dr 



1/2 



— Qk 



3r s 



rs 



(32) 



with k — y at nsco = 3rs = 6GM/c 2 . 

For rin 2> nsco the interface modes are the same as for the Newtonian case. However, as r- m — ► risco the suppression 
effect of cohort [ see Eq. I|26)l] is reduced as the vorticity goes to zero at nsco, so that for r ln = nsco there is no cutoff sound 
speed (see Figures 2-3) for interface mode instability. This is illustrated in Figure 4. Thus if the magnetosphere boundary is 
at nsco, the interface modes will always be present and highly unstable for any sound speed. 



3.4 P-modes with Magnetosphere Boundary 

The boundary condition given by equation (1221! also provides an inner reflection boundary for disc p-modes, which were 
studied in detail in Lai & Tsang (2008). These modes have wavefunctions primarily "trapped" in the wave region between the 
disc boundary n n and the inner Lindblad resonance radius, tilr, where ui — mQ, = —K- Figure 5 depicts an example of the 
p-mode wave function for the same disc model as in Figure 1. The growth rates of these p-modes are determined primarily 
by the outgoing flux at the outer boundary and the effect of the corotation resonance, as discussed in Lai & Tsang (2008). In 
Figure 6 the eigenfrequencies are shown for p-modes in a disc with the density profile E oc r~ p , where p — 3/2 so that wave 
absorption at the corotation resonance is inactive (since in this case the vortensity K 2 /(2f2E) is constant). For the density 
profile p < 3/2, the corotation absorption tends to damp the mode, while for p > 3/2 the corotation absorption enhances it. 
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Figure 5. The p-mode wavefunctions for m = 4, c s = O.lrU, S_/S + = 1/99, = 1 with ui/n in = 1.914 + 0.000127i. The 

real components are shown in solid lines, while the imaginary components are dashed lines. For p-modes the inner and outer Lindblad 
resonances and the corotation resonance are outside rj n , denoted by the dotted lines. 



4 INTERFACE MODES AT THE STAR-DISC BOUNDARY 
4.1 Star-Disc Boundary Condition 

In the case of accretion on to a non-magnetic star, our model consists of a dense uniformly rotating compressible stellar 
atmosphere truncating the accretion disc. This model ignores the structure of the boundary layer. However the qualitative 
properties of the dynamics should be captured for modes with characteristic radial length scale much greater than the radial 
scale of the boundary layer. 

Several studies of CVs (e.g. by examining the rotationally broadened line emissions from the stellar surface) have shown 
that the stellar rotation rates are significantly below the breakup rotation rate (see Warner, 2004), and we limit our examina- 
tions to systems with |f2_ | ^ 0.5f2(ri n +). As in the case of the disc, we consider only the effect of perturbations on a cylindrical 
equatorial surface of the stellar atmosphere (i.e., we are considering a "cylindrical" star). In this region equation (|18[) also 
describes the enthalpy perturbations within the stellar atmosphere. Inside the atmosphere, we assume a small constant density 
scale height, — E/E' = Hs <C r. Equation (I18|l then becomes 



8h" - -^Sti 



2mn D\ c 



0. 



For H-z <C r and Hy, >C c s /Q, this has the solution 
5h oc exp[(r — r in )/H s ]. 

The Lagrangian pressure perturbation at the stellar surface is then 

i8u r 



(33) 



(34) 



AP- = E- 



-'< >/ ■; 



— r(il k - S2 

lOU r 



r(n 2 k - n 2 



i 



i&U r 



(35) 



Once again matching the Lagrangian displacement and pressure perturbation at the interface gives the boundary condition 
for the interface modes for the star-disc boundary case: 



ujSh pc 2 s 


= E_ 


iSu r r 





2 ~ 2 

— +r(Q, -il k ) 



2P.n 

rui 



He 



(36) 



© 0000 RAS, MNRAS 000, 000-000 



Interface Modes in Accretion Discs 9 




0.1 0.2 

sound speed c s 



0.3 



Figure 6. Real and imaginary frequencies for p-modes for the magnetosphere boundary condition with S_ = (l/99)£+ and H_ = 0.5Sli n . 
For £_ <C £+ the p-mode frequencies have very little dependence on f!_ . Here the disc surface density profile is chosen to be S oc r -3 / 2 , 
so that the corotation absorption is inactive and mode growth is purely due to propagation outward at r ou t- 



4.2 Numerical Results 

We repeat the numerical procedure of Section 3 using the radiative outer boundary condition [equation (|24|l ] and using 
equation (|36|l as the inner disc boundary condition. A sample wavefunction for the star-disc interface mode is shown in 
Figure 7. For typical mode frequencies the region just outside the boundary is an evanescent zone; wave propagation becomes 
possible only beyond the outer Lindblad resonance (7"ol)- Figure 8 shows the eigenfrequencies for the lowest order modes with 
m = 1, 2, . . . , 5 as a function of disc sound speed, for representative parameters Hy, = 0.01ri n , E_ = 10E + and f2_ = 0.1f2 + . 
Figure 9 shows the dependence of the mode eigenfrequencies on the density (S_), rotation rate (fi_), and scale height (Ht.) 
of the star. We see that both the real mode frequency uj r and the growth rate Ui do not depend strongly on these parameters. 



4.3 Discussion of Numerical Results 

When the disc is truncated by the star's surface, the effective gravity acts to stabilize the perturbations (since £_ > £+), 
as does the vorticity. Thus compared to the interface mode in the magnetosphere-disc case (Sections 3.1 - 3.2), the mode 
growth rates here are much smaller and are primarily driven by wave propagation through the corotation, beyond the outer 
Lindblad resonance. In the left panels of Figure 8 the eigenfrequencies are shown for typical parameters [E_ = 10E+(n n ), 
0_ = 0.1£2(ri n +), Hy, = O.Olrjn], and disc density index p = 3/2, so that the corotation absorption plays no role. For other 
density indices, wave absorption at the corotation can act to either damp or grow the interface modes (Tsang & Lai 2008; Lai 
& Tsang 2008). For example, in the Shakura-Sunyaev a-disc model the disc solution for the outer disc solution (with free- free 
opacity and gas pressure dominating) has the surface density E tx t- _3//4 , hence the modes would be stabilized by absorption 
at the corotation resonance. However, for models where the disc has density index p > 3/2 at corotation, the corotational 
absorption acts to enhance mode growth, as shown in the right panels of Figure 8. 

The mode eigenfrequencies have very little dependence on the properties of the stellar atmosphere (E_ , Q- , Hy,), as 
shown in Figure 9. The mode frequencies instead primarily depend on the disc sound speed, which in turn depends on the 
accretion rate. Observations of CVs indicate that DNOs are usually only detected in high M states, with the oscillation period 
decreasing with increasing luminosity (Warner 2004). The downward trend of the u r /mQin as a function of c s in Figure 8 
would appear to contradict the observed period-luminosity anti-correlation. But note that in our model, the interface mode 
frequency depends on the sound speed at the inner-most disc region and boundary layer, and it will be necessary to model 
the thermodynamic and radiative properties of the boundary layer in order to compare with observation directly. Also, the 
oscillations of the type considered here would yield periods shorter than the surface Keplerian period, except for the m — 1 
mode. Though higher-m modes would be more difficult to observe due to the averaging out of the luminosity variation, most 
observed DNOs, even those with 1:2:3 harmonic structure (Warner & Woudt, 2005) occur with period roughly at or greater 
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Figure 7. An example wavefunction for a disc/non- magnetized star interface with m = 2, S_ = 10£+, f2_ = 0.1Oj n and Hy, = 0.01rj n 
with the eigcnfrcquency oj/f!i n = 1.378 + 0.0030i. The left side (r < rj n ) of the plot denotes the perturbation eigenfunctions inside the 
star, while the right side shows the disc perturbations. The vertical dotted lines denote the corotation resonance radius (r c ) and the 
outer Lindblad resonance radius (tol)- 




Figure 8. The eigenfrequencies of the interface modes for the star-disc boundary for Hy, = 0.01ri n , S_ = 10E_|_ and f2_ = 0.1Oi n , as a 
function of disc sound speed, for m = 1 . . . 5. The vertically integrated surface density of the disc is given by E+ tx r~ v with p = 3/2 for 
the left panels, so that corotation absorption plays no role, and p = 2 for the right panels, where corotation absorption acts to enhance 
the mode growth. 
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Figure 9. Eigcnfrcqucncies for disc-star interface modes as a function of star density (£-), rotation rate and characteristic scale 

height (£_/£'_ = -He). The canonical values for various parameters arc: c s = 0.2ri n fii n , H_ = 0.1S!i n , = 0.01r-j n and £_ = 10S+. 



than the corresponding surface Keplerian period. These long-period oscillations cannot be explained by the model considered 
here. 



5 CONCLUSIONS 

We have studied the non-radial oscillation modes at the interface between an accretion disc and a magnetosphere or stellar 
surface. Although the models explored in this paper are perhaps too simplified compared to realistic situations, they offer 
some insight into the behavior of the interface modes in various astrophysical contexts. 

Our study of the interface modes at the magnetosphere-disc boundary extended the work by Li & Narayan (2004), who 
considered incompressible disc flow (and therefore could not treat real discs). The model can have very strongly unstable 
modes due to Ray leigh- Taylor and Kelvin-Helmholtz instabilities. In systems where the magnetosphere has developed from 
advection of frozen magnetic flux, the magnetosphere is expected to be roughly rotating with at the Keplerian rate. Since there 
is no shear at the interface, only the Rayleigh- Taylor instability may occur. However the disc vorticty (due to differential 
rotation) acts to suppress the instability, leading to a cutoff below a critical sound speed. Thus a sufficiently hot disc is 
required to generate unstable low-m modes. For magnetospheres rotating with the central star, shearing is expected between 
the magnetosphere and disc, and the Kelvin-Helmholtz instability becomes active. This can help to drive the instability for 
low-m modes to overcome the vorticity in low sound-speed discs. In discs that terminate near the ISCO in a general relativistic 
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potential, the vorticity approaches zero at the inner disc radius, and unstable interface modes can be found for any sound 
speed. 

We can expect a strong dependence of the interface mode growth rates on the sound speed, and thus accretion rate, while 
the real frequencies of these oscillations remain close to (1 — 1.3)mf2i n , and do not depend strongly on the sound speed. It 
is worth noting that the same boundary condition that gives rise to the interface modes also gives rise to intertial-acoustic 
modes (or p-modes) in the disc. 

Although higher-m interface modes are more unstable in our model, these are less likely to be observed due to the 
averaging out of the luminosity variation over the observable emitting region. In addition, if the effect of viscous damping is 
considered (e.g., Wang & Robertson 1985), small wavelength or high-m, perturbations are suppressed. 

These results are for perturbations with no vertical structure, and are applicable mainly to the midplane of accretion 
discs interacting with a magnetosphere. Global 3D numerical studies of Ray leigh- Taylor instability induced accretion onto 
magnetized stars have been performed by Romanova, Kulkarni & Lovelace (2008) and Kulkarni & Romanova (2008), and show 
such small m instabilities in the disk midplane. The low-m oscillations at the magnetosphere-disc interface may be relevent 
to the high-frequency QPOs observed in some NS and BH X-ray binary systems (Li & Narayan, 2004; see Section 1 of Lai & 
Tsang for a critical review of various theoretical models), although to obtain the correct QPO frequencies for the BH systems, 
the disc inner radius must lie outside the inner-most stable circular orbit. 

For the star-disc boundary model considered in in this paper, the interface mode growth rates are much smaller than for 
the magnetospheric case, since the effective gravity now acts to stabilize the system, and the Rayleigh- Taylor instability is 
inactive. The modes discussed here are unstable due primarily to propagation through the corotation. With sufficiently steep 
disc density profile (E oc r~ v with p > 3/2), corotation absorption can also help to drive these modes, as studed previously by 
Tsang & Lai (2008) and Lai & Tsang (2008). Such modes may be responsible for the high-frequency (of order the Keplerian 
frequency at the stellar surface) dwarf nova oscillations observed in CVs, although oscillations with longer periods would 
require a different explanation. 
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APPENDIX: PLANE PARALLEL FLOW WITH A COMPRESSIBLE UPPER LAYER 



Consider a system consisting of two fluids in the gravitational field g = — gz. The upper fluid (z > 0) has density p = p+e~ Z//fl * 
(with H z — c 2 /g, where c s is the sound speed) and horizontal velocity u+ along the x-axis; the lower fluid (z < 0) is 
incompressible with density p_ and horizontal velocity U-. 
The linear perturbation equations for the upper fluid are 



jW- (*» + »*,) = 
at 

— Su + (u- V)<5u+ (<Su- V)u D = -VSh 

where Sh = 5P/ p. For perturbations of the form e 1 
d 

— iiuSp + ikpSu x + —(p6u z ) = 0, 
oz 



these become 



i(Z>8u z = 
iuj8u x = 







Sh, 



dz 
ikSh, 



where u)+ = uj — ku+. Assuming the perturbation is isothermal, so that 5P = c^Sp, we obtain 
1 



5h"(z) 



■5h'{z) - (k 2 - u + /c 3 )8h(z) = 0. 



The two independent solutions of equation ((42 



Sh oc exp 



z 
2Hl 



1 ± sf\ + AHlkl 



where k\ = (k 2 — Hi 2 /(?). Obviously, the physically relevant solution is 
Sh oc exp(— kz), 



(37) 
(38) 

(39) 

(40) 
(41) 

(42) 

(43) 
(44) 
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with k — (VI + 4H|fc| — 1)/2H Z . This gives, by equation (|40|1 the Eulerian pressure perturbation for the upper region: 
SP+ = -!^± Suz+ . (45) 

For the lower region the fluid is an incompressible potential flow with 8u — V6ip with Sip satisifying \7 2 8ip — 0. For z< 0, 
the appropriate solution is 

dtp oc exp(fcz). (46) 

This gives V[iu}5ip + 5P_/p_) = 0. The Eulerian pressure perturbation in the lower region (z < 0) is then: 

5p= ip^- 5uz _ (47) 
k 

Matching the Lagrangian displacement and Lagrangian pressure perturbation across the boundary between the upper 
and lower regions we get 

~ 2 ~ 2 

-±f±+p + g = - —+p-g, 48 
k k 

which can be written as the quadratic: 



p+~2ku+ + p-2ku- 



P+T + P- 



p+jk 2 u\ + p~k 2 u 2 _ p,-p 

k : t , + g „ P+ fc P : ■ (49) 



P+T + P- P+I+P- 



where k has a non trivial dependence on uo. This has a solution 



_ = k(p+u+ + p-u-) ± / k 2 (u+ - u-)p+p- kg(p+-p-) 

p++p- V (p++p-) 2 p++p- 

where p + — p + k/k. 
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